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Abstract 

We critically examine recently published results on the thermal 
production of massive vector bosons in a quark-gluon plasma. We 
claim the production rate is a collinear safe observable. 
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*and a recent unpublished paper by Wong. 



In two recent papers Kapusta and Wong [1,2] calculate the order g1 
finite temperature corrections to dilepton and Z boson production in a hot 
quark-gluon plasma in the limit where the dilepton or Z boson mass M is 
large compared to the temperature T. This involves calculating the two-loop 
QCD corrections to the imaginary part of the virtual photon or Z self-energy 
diagram. The result, in both cases, is found to be 
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where k^ is a "cut-off on the four momentum transfer carried by the ex- 
changed quarks in the plasma" . This cut-off is identified with "the effective 
mass of a quark propagating through the plasma with a typical thermal 
momentum" {m^^ = g^T'^/S). If ttIcs is taken to vanish then the two-loop 
expression diverges. This signals a collinear singularity shielded by the non- 
vanishing thermal mass of the quark. 

The result eq. (1) above is very surprising as it contradicts a series of 
works, published over the last twelve years, which have been ignored by 
the authors of refs. [1,2]. Indeed the production rate of a heavy particle 
in a thermal medium has been discussed several times in a variety of the- 
ories (QED, QCD, scalar theories) and has been found to be an infra-red 
as well as a collinear safe observable [3]. These studies have been carried 
out in the real-time formalism using the full 2x2 matrix formalism and the 
imaginary part of the relevant two-point function has been calculated using 
the Kobes-Semenoff cutting rules [4]. All works have the common feature 
that intermediate steps in the calculation require the introduction of regu- 
lators which are taken to zero after all pieces are put together to construct 
the physical observable. The results, up to three loops in some cases, have 
always been found to be finite after all regulators are taken to zero. 

In refs. [5 — 7] the very same process as in ref. [2] is considered namely, 
lepton pair production in a quark gluon plasma (with massless quarks and 
gluons at zero temperature). In [5] two independent regulators, a gluon 
mass and a quark mass, are introduced; in [6], instead of introducing mass 
regulators, the calculation is carried out in dimensional regularisation (n = 
4 — 2e) ; in [7] the calculation is done in four dimensions but a small mass is 
given to the gluon. In all these cases it was found that the imaginary part 
of the photon two-point function, evaluated in the two-loop approximation, 
is finite when the regulators vanish, indicating the absence of infra-red and 



collinear singularities for this observable. In [8], the same calculation is 
repeated in n dimensions but keeping the quark massive: no logarithmic 
sensitivity to the quark mass is found in the final result. It is interesting to 
note that Cleymans and Dadic [9] have calculated the imaginary part of the 
photon two-point function for a space-like photon (M^ < 0) and they also 
found a finite result when the regulators are vanishing. Finally the problem 
of infra-red and collinear singularities was examined in a scalar theory [gcj)^ 
in six dimensions). Cancellation of all singularities was found in the two- 
loop approximation [10] as well as in the three-loop approximation [11, 12]. 
All these results are therefore in contradiction to eqs. (1) (see also [13, 14] 
for work on related processes). 

More work has been performed concerning the production of a virtual 
photon in a quark- gluon plasma and the 0{g1) terms have in fact already 
been explicitely calculated ten years ago! Independently of the regulari- 
sation method used (dimensional regularisation in [15] and gluon mass in 
[7]) the 0{g1) corrections have been found to be identical as it was checked 
in numerical calculations [7,15]. Furthermore, in [15] a simple algebraic 
expression is derived in the limit of interest in refs. [1, 2] (M ^ T) 
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(see eq. (4.5) in [15]) where the factor 1 in the right-hand side is just the 
0{g1) correction factor for e+e~ annihilation into hadrons at zero temper- 
ature^ and all terms in T"^ /M"^ vanish leaving terms in T'^/M^ as leading 
thermal corrections. On the other hand, a logarithmic behavior appears in 
the hmit M <^T, 

It is interesting to note that the above equation has been obtained in an 
explicit calculation in the bare theory [15] as well as in the limit of the 
effective theory [16] when the intrinsic thermal regulators (quark and gluon 
thermal masses) are becoming vanishingly small [17]. 

Not enough details are given in refs. [1,2] to locate the origin of the 
error. In a recent preprint however, Wong [18] presents the formalism upon 
which eqs. (1) are based and he suggests that the cause of the discrepancy 

'^This factor 1 can be seen as a check on the calculation and arises from the process 
qqG -^ 7* after summing over all allowed phase space of the initial quark, anti-quark and 
gluon. 



between refs. [1, 2] and previous works is due to the fact that previous works 
("for example [15]") include a form of resummation. This is incorrect since 
all previous works [5 — 15] are based on a strictly perturbative approach 
similarly to the approach followed in [1,2]. 

Starting from the imaginary-time formalism the author of [18] looks for 
an expression of the imaginary part of the photon two-point function. He 
decomposes it into two pieces: a "better known" I{k) part and a "not well- 
known" J{k) part. A physical interpretation is proposed for each of these 
pieces. The I{k) part is obtained "by putting the three internal lines in 
the two-loop graphs on shell" and corresponds to physical processes evalu- 
ated in the tree approximation. They are labeled as Compton scattering, 
decay with photon emission and vector-photon fusion and they are the pro- 
cesses which resemble those obtained "by cutting rules at zero temperature" . 
They appear with the appropriate thermal Bose-Einstein {f^^' or 1 -|- f^^') 
or Fermi-Dirac (/'"-' or 1 — /'"-*) distributions depending on whether the 
particles are absorbed in the initial state or produced in the final state. The 
J{k) part is obtained by putting two of the three internal lines on-shell while 
the uncut internal line is attributed, rather arbitrarily, the thermal weight 
1/2 -|- /*•+) or 1/2 — /("' for a boson or a fermion respectively. These terms 
are labeled "interference terms". Let us remark that the thermal weights 
appearing in I{k) and J{k) should automatically come out from the ther- 
mal cutting rules and should not be the result of an educated guess. Strictly 
adhering to the cutting rules is the approach followed in refs. [5 — 15] and 
this leads to the correct result. There, the I{k) type terms and the J{k) 
type terms were both included with the appropriate thermal factors and 
they were labeled "real diagrams" and "virtual diagrams", respectively, by 
analogy with the zero temperature case. 

Perhaps the best way to gauge the approach followed in [18] is to compare 
it with the retarded/advanced version of the real-time formalism [19 — 22]. 
Indeed this formalism allows the construction, in the real-time approach, 
of Green's functions which are directly comparable to those obtained in 
[1,2] from the analytical continuation of the imaginary-time Green's func- 
tions. Expressions for one and two-loop Green's functions are explicitely 
constructed [19,23]^. More precisely, the very same expressions needed for 
the calculation of [1, 2] can be found in [24] in the slightly more complicated 
case of the resummed theory (see eqs. (24) and (29) of [24]). Going back to 
the bare theory and comparing to the expressions in [18] should be an easy 

^A brief discussion is also given on the choice of the thermal factors of the form (a±/' ') 
where a is an arbitrary constant. 



task. 

Another subtlety in thermal calculations, which has in the past led to dis- 
crepancies between different calculations, is related to the treatement of the 
fermion effective mass in a hot medium. Thermal corrections to a fermion 
propagator do not break chiral invariance [25 — 27], i.e. the thermally gen- 
erated mass arises from a term of the form ^/^ where I^ is a vector-like 
function depending on the momentum of fermion ^ and on the temperature. 
In actual calculations this leads to distinguish [28 — 31, 15] between a "kine- 
matical mass" which incorporates the thermal effects and which appears in 
the phase-space constraints and a "dynamical mass", the fermion mass at 
zero temperature, which occurs in the evaluation of the matrix elements (i.e. 
the thermal mass shift drops out in the trace evaluation). Using a scalar 
mass shift 6m^. ,^^ in the calculation would lead to incorrect estimates 

i ncrmal 

of the production rate [30] . 
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